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Abstract: With the probe limit, we investigate the behavior of the electric permit- 
tivity and effective magnetic permeability and related optical properties in the s-wave 
holographic superconductors. In particular, our result shows that unlike the strong cou- 
pled systems which admit a gravity dual of charged black holes in the bulk, the electric 
permittivity and effective magnetic permeability are unable to conspire to bring about 
the negative Depine-Lakhtakia index at low frequencies, which implies that the nega- 
tive phase velocity does not appear in the holographic superconductors under such a 
situation. 



Contents 



|1]. Introduction 

0. Holographic model of superconductors 



2.1 Abelian-Higgs model and holographic dictionary 



2.2 Holographic normal phase and superconducting phase 



Refractive index in holographic superconductors |5| 
5TT] Holographic setup [5| 



3.2 Numerical results 



Conclusions |12 



1. Introduction 

Since the advent of AdS/CFT correspondence, not only do the main efforts include 
the examination of its validity case by case, but also involve its applications in various 
fields, in which condensed matter physics has recently received a special attention. The 
motives are twofold. On the one hand, there exists many strong coupled systems in 
condensed matter physics, which are intractable by the traditional approaches. While 
AdS/CFT correspondence, as kind of strong/weak duality, can provide a powerful tool 
to attack them. On the other hand, unlike other disciplines such as particle physics and 
cosmology, condensed matter physics allows one to engineer matter such that various 
vacuum states and phases can be created in the laboratory, which may in turn provide 
the first experimental evidence for AdS/CFT correspondence. For a review of this 
subject, please refer to [Jj, 0, [|, |J. 

In particular, most recently such a bulk/boundary duality has been applied to 
explore the optical properties of strong coupled media which admit charged black holes 
as a dual gravitational description and it is shown that this type of media generally 
have a negative refractive index at low frequencies 0. Along this line, this paper is 
a first attempt to investigate how the refractive index and related optical properties 
behave in the holographic superconductors and to see whether the negative refractive 
index shows up at low frequencies. 
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The gravity dual of superconductors was firstly established to model the s-wave 
superconductors 0, [F). Later both the p-wave and d-wave holographic superconductors 
were also realized [||, ^ |10 , [TTJ, [12], |T3|. Please refer to |14] and references therein for a 
review of various properties related to such holographic superconductors. 

In the next section, we shall recall the holographic model of s-wave superconductors, 
where to make our life easier we would like to work in the probe limit and focus on the 
special case in which the complex scalar field is massless. After the setup, in Section || 
we will present our numerical results for the behavior of the relevant optical quantities 
in such holographic superconductors. Conclusions and discussions will be addressed in 
the end. 



2. Holographic model of superconductors 
2.1 Abelian-Higgs model and holographic dictionary 

The bulk dual contents of a holographic superconductor with a s-wave order parameter 
in 4+1 dimensions include the gravitational field, the U(l) gauge field, and the complex 
scalar field of mass m and charge q. The corresponding bulk action is given by[|TT|j 



'bulk 



where F = dA, and D = V — %A where V is the covariant derivative compatible with 
the metric. In what follows, we will work in the probe limit, i.e., ^ — > in which 
gravity decouples from the matter fields and the backreaction from the matter fields 
can be ignored. Thus a solution to Einstein equation is a black brane, i.e., 

7-2 j 2 

ds 2 = —[-f( u )dt 2 + dx 2 + dy 2 + dz 2 + -^-], (2.2) 
u 1 f(u) 

where f(u) = 1 — (^-) 4 with Uh > the horizon. Obviously, in the limit u — > 0, such 
a solution asymptotically becomes anti-de Sitter. Next around this background the 
equations of motion for the matter field are given by 

= D a D a $ - m 2 $ = -;L(0 a - iA a )[^g a \d b - iA b )<&\ - m 2 $, 



-g 

= V a F ab - i($D b <$> - $Z) 6 <i) = -^=d a {^F ah ) - i[${d b - iA b )$ - $(d b + iA b )$]. 



(2.3) 



By the holographic dictionary, the boundary value of the metric acts as the source 
for the energy momentum stress tensor T^ u on the boundary, the bulk gauge field A^ 
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evaluated at the boundary serves as the source for a conserved current associated 
with a global U(l) symmetry, and the near boundary data of the scalar field $ sources a 
scalar operator O with the conformal scaling dimension A satisfying A (A — 4) = m 2 L 2 . 
In what follows, we shall be focusing on the particular case, i.e., m 2 = 0, which yields 
an operator of dimension four by normalizability. Whence the asymptotic solution of 
$ and can be expanded as 

$ = _)L($(°) + $(V + •••), 

A, = Af + i[^V - g^d^Af - d,Af)u 2 In ^ + • • • (2.4) 

with e< 1 a UV cutoff and = lim^o j^g^v Then it follows from the holographic 
dictionary that the expectation value of the corresponding boundary quantum field 
theory operators (O) and ( J M ) can be obtained by variations of the renormalized action 
with respect to the sources, i.e., 

(O) = lim 1_ 5 4- =$W 

(J") = lim — L=j^ = g^ v [AV + cg^d^Af - ft,^)], (2.5) 

where the constant c is renormalization scheme dependent [jDj] , to be fixed in the later 
discussions. 

2.2 Holographic normal phase and superconducting phase 

With the following ansatz, i.e., 

$ = A* = 0(z), A x = 0, A y = 0, A 2 = 0, (2.6) 

the equations of motion can be reduced to 

f 3 (h 2 
= $'' + (!- 

/ u f 2 

1 07-2(1)2 

= 0"- -<t>' - ^-H, (2.7) 

where the prime denotes the differentiation with respect to u and $ has been taken 
to be real by taking into account the fact that the u component of Maxwell equations 
implies that the phase of $ is independent of u. Note that the regularization condition 
requires <p = on the horizon. Then multiplying the first equation in (|2.7|) by /, one 
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Figure 1: The condensate as a function of charge density, where the temperature is fixed to 
be 1 as the reference scale. 

can find $' = on the horizon. Thus for the above two second differential equations, 
we have only a two parameter family of solutions, which can be labeled by $>(uh) and 
4>'(uh)- For convenience, we will set q — 1, L — 1, — 1 in the later calculations. The 
result for any other can be easily restored by the following scaling rules, i.e., 



$ (0) K) 




$ (1) K) 






< 


4 0) K) 


= -4 0) (i), 


4 1} K) 




T(u h ) 


1 



U h 7T 



Firstly, it is easy to find the trivial solution as 

$ = O,0 = C(1-m 2 ), (2.9) 

which corresponds to the normal phase on the boundary with the chemical potential C. 
On the other hand, according to the numerical calculations, by fixing the charge density 
or chemical potential the non-trivial solution emerges below a critical temperature 
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Figure 2: The condensate as a function of chemical potential, where the temperature is fixed 
to be 1 as the reference scale. 



T c , and the condensate increases when the temperate is lowered [16] . Its boundary 
correspondence is the superconducting phase. 

By the above scaling rules, such a non-trivial solution can also be thought of as 
the superconducting phase emerging above a critical charge density p c or a chemical 
potential S c if the temperature is fixed. We would like to plot the corresponding 
condensate as a function of charge density and chemical potential in Fig.|l] and Fig.|| 
respectively. 



3. Refractive index in holographic superconductors 



3.1 Holographic setup 

Before proceeding, we would like to introduce some quantities we want to calculate and 
their relations without going into details. In the linear response theory, the induced 
electromagnetic current J is related to the external vector potential A as Ji(u,k) = 
Gij(oj,k)Aj(uj,k), where G is the retarded Green function. On the other hand, in 
the Laudau-Lifshitz approach to electrodynamics of continuous media[[L7|, [18], [19|], for 
the isotropic media, the transverse part of the dielectric tensor is determined by the 
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transverse retarded Green function as follows 1 1| |20] 



47T 

e T (u,k) = 1 + — G T (u,k), (3.1) 

from which the corresponding electric permittivity and effective magnetic permeability 
can be expressed as 

An 

e(u) = 1 + —G T0 (u), 
ur 

= i-Ja^y (3 ' 2) 

Here Gtq{oS) and Gt2{^) are the expansion coefficients of the transverse retarded Green 
function in k, i.e., 

G T (u, k) = G T0 (u) + k 2 G T2 (u) + ■■■ (3.3) 

with Gt{u), k) = (5^ — K^-)Gij(uj, k). With this, the refractive index can be given by 
n 2 (cu) = e(u)fi(u). But to identify whether our holographic superconductors display 
the opposite phase velocity to the power flow, we shall appeal to the Depine-Lakhtakia 
index, which is defined as 

n DL {u) = \e(u)\Re\p,(u)] + \fj,(u)\Re[e(co)}. (3.4) 

As shown in |Tl , the phase velocity is opposite to the power flow if and only if Udl < 0. 

Now let us move on to the strategy to calculate these quantities by holography. 
By symmetry, in order to consider the refractive index in the above holographic su- 
perconductors, it is enough to consider the transverse fluctuations of the field A as 
follows 

SA X = a x (u)e~ My . (3.5) 

Note that such a fluctuation decouples and has no back reaction, satisfying the equation 
of motion 

/ » ' P f M 
Whence a x goes like (1— u) ±l i near the horizon, corresponding to the outgoing and ingo- 
ing boundary conditions respectively. We here choose the ingoing boundary condition, 
which will lead to the retarded Green function in the dual field theory ^2|. Speaking 
specifically, with such an ingoing boundary condition, we assume the corresponding 
asymptotic expansion a x goes like 



= < + (V - ~K + (—-_- —)a x . (3.6) 



a(°) + l[4 1 V-(a; 2 -/c 2 )4 Vln^ + --. (3.7) 



1 Note that the expression in (E) is different from ours due to the different conventions used for the 
retarded Green function. 
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where we have used Eg . ( |3~5l) and the second eguation in (|3.7|) . By Eq. fl2.5Q , the induced 
current is given by 

j x = + c(u 2 - k 2 )a x °l (3.8) 
Whence the retarded transverse Green function can be obtained as 

G T (u,k) = — +c (u 2 -k 2 ). (3.9) 

An analytic solution of Eq.( [3.6|) does not appear to be available. So in the subse- 
guent section, we shall resort to Mathematica for numerical calculation, where in order 
to obtain Gto(oj) and Gt2(w), it is advisable to adopt an alternative approach by firstly 
expanding a x in series of k, i.e., 

a x = a x0 + k 2 a x2 H , (3.10) 

from which Eq.( [3.6|) becomes 



= a x0 + -7 - - Ko + (72 ~ Oa;0 ' 
f u f 2 fu 2 

,f' ^ / V 2$2 N 1 

= a x2 + (— - - )a x2 + (— - -7-3)0x2 - 7^0- (3.11 



Near the boundary, the asymptotic expansion for a x0 and a x2 can be read out of Eg. Q3.7| ) 

as 

a^o - Oxo + 2^ a xo u ~ u a xo u ln -J H 

(0) , If (1) 2 1 2 (0) (Ok 2 1 W i , /, i n \ 

«x2 = <2 + 2^2" - (w < 2 ; - < > In -j + • • (3.12) 

Whence we have 

(1) (i) (0) 

G T^) = a -W ) [^f ) - a -^]-c, (3.13) 

a x0 a x0 a x0 

where c can be fixed by reguiring that Gto approaches zero at large freguencies as the 
system has no time to respond to the rapid variation of external fields. 
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Figure 3: The imaginary part of the ey for real uj and k at temperature ijr = 0.85. 
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Figure 4: The imaginary part of the €t for real uj and k at temperature ^- = 0.25. 
3.2 Numerical results 

Here we plot all the results by fixing the charge density to be 1 but varying the tem- 
perature with respect to the critical temperature and consequently the condensate. 
As explained in || , firstly we have to ensure that the system is in thermodynamical 
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Figure 5: The real part of the electric permittivity as a function of frequency with various 
temperatures. 

Im[e] 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 T c 



Figure 6: The imaginary part of the electric permittivity as a function of frequency with 
various temperatures. 

equilibrium, i.e., the imaginary part of Gt, or equivalently the imaginary part of ej- 
should be greater than zero for real u and k. This is actually guaranteed by the holo- 
graphic duality, as the bulk background is static with the black brane. For illustration, 
we would like to plot the imaginary part of in Fig.|] and Fig.| to demonstrate that 
this is the case. 

Then we plot the real and imaginary parts of the electric permittivity in Fig.|| 
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Figure 7: The real part of the effective magnetic permeability as a function of frequency 
with various temperatures. 



and Fig.[5] respectively. By the scaling rules, to lower the temperature is equivalent to 
increase the charge density and chemical potential. Therefore the electric permittivity 
here demonstrates the similar behavior at low frequencies as that considered in ||. In 
particular, the larger are the charge density and chemical potential, the larger is the 
regime for the negative real part of the electric permittivity. In addition, below the 
critical temperature, both of the real and imaginary parts of the electric permittivity 
diverge at u — 0. By the fact that the electric permittivity is related to the conductivity 
as e(cj) = 1 + i^j-a(oj), the real part of the electric permittivity diverges as — -tt, and 
accordingly the imaginary part blows up like due to the Kramers- Kronig relations 
from causality . 

Next we plot the real and imaginary parts of the effective magnetic permeability in 
Fig.[F] and Fig.^] respectively. The real part of the effective magnetic permeability still 
displays the similar behavior at low frequencies as that in || , while the imaginary part 
shows a totally different behavior, i.e., there is no blow-up at uj = 0. 

Such a difference eventually leads to the conclusion that at low frequencies the 
negative phase velocity does not show up in the holographic superconductors at least 
within the accuracy of our numerics, which is illustrated by the non-negative Depine- 
Lakhtakia index in Fig.§. In hindsight, the reason may arise in the fact that the bulk 
background is essentially kind of neutral black hole although the electromagnetic field 
has an influence on the fluctuation equation Q3.6|) indirectly through the condensate. 

At last, we would like to plot the real and imaginary parts of the refractive index 
in Fig.ffO and Fig.O respectively. Note that the ratio of the imaginary part of the 
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Figure 8: The imaginary part of the effective magnetic permeability as a function of fre- 
quency with various temperatures. 
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Figure 9: The Depine-Lakhtakia index as a function of frequency with various temperatures. 



refractive index to its real part describes the ratio between dissipation and propagation. 
We see that the there exists a characteristic frequency, below which light can not 
propagate in our holographic superconductors while above which light can propagate 
without dissipation. 
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Figure 10: The real part of the refractive index as a function of frequency with various 
temperatures. 
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Figure 11: The imaginary part of the refractive as a function of frequency with various 
temperatures. 

4. Conclusions 

We have carried out the numerical analysis of the optical properties of the s-wave 
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holographic superconductors in the probe limit. In particular, we have calculated 
the electric permittivity, effective magnetic permeability, refractive index, and Depine- 
Lakhtakia index by the holographic duality, which is generically difficult to be achieved 
by other approaches. As a result, the electric permittivity and effective magnetic per- 
meability conspire to make our holographic superconductors robust against a negative 
phase velocity at low frequencies. 

We conclude with some generalizations of our work in various directions. Firstly, 
although the cases for other scalar masses are expected to demonstrate the qualitatively 
similar properties as the massless case, it is interesting to investigate how the specific 
optical properties of holographic superconductors depend quantitatively on the mass 
m. Secondly, it is natural to extend our work to the 2+1 dimensional holographic 
superconductors^]. More importantly, as the superconductors are believed to have 
the potential to support low losses, which is critical for many applications like super- 
resolution imaging, it is rewarding to search for the negative phase velocity in other 
situations. In particular, it is worthwhile to go beyond the probe limit to explore 
the full behavior of the theory, where the negative refractive index is expected to 
appear at low frequencies because the background is essentially kind of charged black 
hole. In addition, as suggested by the early experimental implementation of negative 
refractive index in the ferromagnet-superconductors ||24j| , it is highly possible to see the 
negative refractive show up when the external magnetic field is added to our holographic 
superconductors. Furthermore, it is also intriguing to extend our analysis to both the 
p-wave and d-wave holographic superconductors, where some new features should come 
out due to the nonisotropy of the dual media. Finally, it is also tempted to generalize our 
work to some non-minimal models of holographic superconductors [p5|, [26|, E7I, EM |29|, RU 
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